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We investigate the conditions for entanglement in a system of two atoms and two photon modes in 
a vacuum, using the Jaynes-Cummings model in the rotating-wave approximation. It is found that 
the strength of atom entanglement is a periodic function of time, generalizing the results of other 
workers. We explicitly show that our results are in agreement with existing results and reproduce 
existing entanglement conditions under appropriate limits. Results for the two-atoms two-photons 
system are generalized to the case of arbitrary values for the atomic energies, relative to photon 
frequencies. We consider both the resonant and non-resonant conditions. 

PACS numbers: 



I. INTRODUCTION 



Entanglement of quantum states is not a new concept; however it is not a property of Fock Space [1]. Therefore it 
does not appear automatically in a vacuum and one has to develop a special representation using second quantization 
to entangle atomic states with vacuum. 

This phenomenon is still not well-understood. There are only a few existing references [1-4] in the literature that 
have studied the possibility of entanglement in the second quantization formalism using simple theoretical models. 
Pawlowski and Czachor (PC) in Rcf. [1] have used a simple model in a system with two atoms and two photon modes. 
They found that the entanglement of two atoms with the vacuum can occur using the canonical commutation relations. 

On the other hand, the Jaynes and Cummings (JC) model [5] is considered to be one of the most appropriate models 
for the purpose of analyzing ion traps in cavities. The JC model, being a nonlinear model, gives a good theoretical 
tool to study ion trapping in a cavity using quantum electrodynamics. 

Hussin and Nieto (HN) in Ref.[6] have studied the JC model in the rotating wave approximation (RWA) to construct 
coherent quantum states using ladder operators. We use the same model in the same rotating wave approximation to 
study entanglement of more than one atom with vacuum. However, we consider a more general case of two-atoms and 
two photon modes. For this purpose we use a form of the JC hamiltonian given by HN and several other references 
therein. The stationary states of the JC model in the rotating wave approximation are given in Ref. [6]. 

The master equation for the cavity losses can be taken as [7] 

dpj i ( 4.I4- lf\ 

-j£ = T 7/ • " ["jI'j", - ^"j"<I>j - 2 p i a i a n ' 

for system j, where j = A or B. The model in the present study assumes two atoms labeled A and B, and also two 
distinct photon mode, also labeled as A and B. Atom j interacts with mode j only, pj is the density matrix of the 
atom-cavity system for the jth mode and jth atom. The factor 7 in the second term is the rate of loss of photons from 
the cavity, due to imperfect reflectivity of the cavity mirrors. Hj is the JC hamiltonian for the jth particle, given as: 

Hj ^fiujjiNj + -)I + -Eja z:j + /iKj{afa-j + a 3 a +j ). (1) 

I is the identity matrix. The Hamiltonian, Eq. (1), is identical to that used in Ref. [7], except that we have included 
the zero-point energy contribution to the photon energy. 

In Eq. (1), each atom, (j = A, B for atoms A and B) has a ground state |— >j and an excited state |+ >j as shown 
in Fig.l. 
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FIG. 1: Absorption of photon in an atom 



The atomic energy difference for atom j is Ej . The frequency of photon mode j is denoted as Wj . The strength of 
the interaction between atom j and mode j is Kj. The x, y and z components of the matrix a are the Pauli matrices. 
The Pauli matrices act on the atomic states exactly as the Pauli matrices act on a spin one-half particle, with the 
lower row of the spinor representing the atoms ground state, and the upper row representing the excited state. We 
also adopt the usual definitions of the raising and lowering operators: <r + = a x +ia y and er_ = <r x —i<j y . The operators 
for the photon mode are defined in the usual way. The operator aj is the destruction operator for photons in mode 
j. Also, Nj = dj^cLj is the number operator for mode j. The choice for the zero of energy of atom j is taken to be 
midway between its ground state and excited state energies. 

Moreover, 



a +l = \+>j< ~\j 

<?zj = 1+ >j< +\j ~ I" >j< ~\j (2) 

The operator a Z j acts on atomic states of atom j as the Pauli spin matrix in the z direction, with the excited state 
considered as the 'up-spin state' and the ground spin state is taken as the 'down-spin state'. Then the second term 
of Eq.(l) denotes the unperturbed atomic state energies. 

The first term of Eq.(l) is the unperturbed photon energy of mode j, with jfuvj as the mode energy. Ej is the 
change of atomic energy due to the absorption of a photon in mode j. 

Our main goal is to study the probability of entanglement if a photon in mode A or B is absorbed or emitted by 
an atom in state |+ > or |— >. These atoms may or may not be identical. We shall specialize to the case of identical 
atoms later, for the sake of simplicity. 

We use a straightforward generalization of the Hamiltonian used by HN[6], extending it to a two-atom, two mode 
system. The complete Hamiltonian is a sum over j (j = A, B) of Eq. (1). 

We present our calculations of the stationary states in Section II. In Section III, we study the time evolution of 
quantum states. Section IV compares our results with those of PC. Finally, section V is devoted to the discussion of 
results and of possible technical applications. 



II. STATIONARY STATES 



We start our calculations with the Hamiltonian of Eq. (1) for the system j, (j = A, B), which represents the JC 
model in the RWA. 

Each photon mode j can cause a transition of atom j between its ground state and its excited state via the emission 
or absorption of a photon. 

It can be shown in a strightforward way that the total number of excitations in the cavity-atom system are given 

by 



This is a constant of motion for the jth mode. From this, one can easily obtain the eigenvalues and eigenstates [7]. 
The Hamiltonian H for the two atoms, two modes problem can be written as 



H = Ha + Hb 



(4) 
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In writing Eq. (4), it is assumed that the two photon modes A and B are distinct. It can be noted that the model 
expressed by Eq. (4) could be extended to include an interaction between atom i and mode j for j =/= i . This would 
be a generalization of the problem studied in references 1 and 2. It is also worth-mentioning here that Eqs.(l) and 
(4) contain the unperturbed atom and mode energies( first two terms of Eq.(l)), as well as the interactions. In the 
study of PC, the Hamiltonian includes only the interaction terms, the last term in our Eq.(l). Following HN [6], we 
introduce the dimensionless parameters A,,- and tj by 



/iu>j =E j (l + e j ) 



(5) 



and 



/iKj 



(6) 



Y 



The parameter €j is a " detuning parameter" in that it is a measure of the deviation of the photon energy /hujj 
from the atomic energy difference Ej. Then Eq.(l) can be rewritten as : 



Hi 



(1 + e^EjiNj + I) + -J^l + XjEjiuja j + a j( j +j ) 



(7) 



The lack of coupling between Hamiltonians Ha and Hb means that a complete set of stationary states of H can be 
formed from products of the stationary states of Ha and Hb ■ The ground state of Hj is simply 



\G>< 



10; 



(8) 



where |0; — >jdenotes the state in which the photon state is the vacuum and atom j is in its ground state. Its energy 
is (Eg), =ej Ej/2. 

The normalized excited states of Hj can be enumerated by n= 0, 1, 2, They are 



l^n >j= (cos9 nj )\n;+ >j -(sm6 nj )\n+l;- >j 



(9) 



and also 



>j= {sme nj )\n;+ >j + (cos 9 nj )\n + 1;- > 3 



(10) 



with energies 



Etj = (1 + + l ) E j ± + 



(11) 



In Eqs. (9) and (10), \n; ± >j denotes the state with atom j in atomic state |± > and with n photons in mode j . 
The angle 6 n j appearing in Eqs. (9) and (10) is defined by 



cos# nj = \j(q n +i,j + y)/ 2 9n+i,j 



(12) 



Finally, q n j is defined by 



sin 6n 



( iv W( " 



n+l,j 



■)/2g n 



+ 1,3 



(13) 



1n, j = \l(^) + nX 2 j 



(14) 
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To write down a basis of stationary states for the full Hamiltonian H, we only need to take products of the stationary 
states of systems A and B. Then, the ground state of H is 



The excited states are 



\G > =\G> A \G> B 
= |0;- > A |0;- > B 

l<A) >a \G >b, 



(15) 



.(a) 



and 



\ G >A Wn(B) >B,- 



■(b) 



\^n(A) >A \i>n{B) >B, 



.(c) 



(16) 



with n(A) and n(B) each taking values 0,1,2, In Eq.(16c), all four choices of the signs + and - must be included. 

The problem studied by PC focuses on the vector space V spanned by the four states 



$1 


> 


= |0; 


- >a |l; 


- >B 




> 




- >a |0; 


- >B 


$ 3 


> 


= 10; 


+ >a |0; 


- >B 


$ 4 


> 


= 10; 


- >a |0; 


+ >B 



(17) 

These four states of Eq.(17) are the tensor product of atomic ground states and excited states with known photon 
modes. These states can ultimately show entanglement. 

The study of PC considers the choice of the initial state (t=0) as 



\1>a >= (^=)(|*1 > +|*2 >) 



(18) 



and the time development of \ip a > is analyzed. PC has studied the resonant case (€= 0)). They have only 
employed the interaction term of the JC hamiltonian (Eq.(l)). We have included the non-resonant case in the next 
section also. 

III. TIME EVOLUTION IN JC MODEL 
A. General Results 

Eqs.(15) and (16) give the stationary states for the system of two atoms and two photon modes. Inspection of this 
equation shows that the vector space V is also spanned by the following four stationary states: 



IV>i > 

\lp2 > 

IV'S > 

|V>4 > 



\G > A \lp >B 
\i/Jq >A \G > B 
\G > A IV'o" >B 
|-0O" >A \G > B 



(19) 



If the initial state is any state in V, its evolution is found by expansion of the initial state in the set \ipk >, where 
k=l, 2, 3, and 4. If each term in the expension is multiplied by exp —iEkt/ /i, (with Ek equal to the energy of \ipk >), 
we have the evolution of the initial state. 

The energies Ek can be obtained from Eq.(ll) by adding the energies of systems A and B. The results are 
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El 


= E GA + E QB = 


£aE a 
2 


+ (1 +€ B )E B 


— <1obE b 


E 2 


= E GB + E~ A = 


£bEb 
2 


+ (1 + e A )E A 


- QoaE a 


E3 


= E GA + E+ B = 


£aE a 
2 


+ (l + e B )E B 


+ QobE b 


£4 


= E GB + E+ A = 


€ bEb 
2 


+ (1 + e A )E A 


+ QoaEa 



(20) 

We now consider the case for which the initial state is given by Eq. (18). The expansion of \ip a > in the states 
\tjj m > is 

! 4 

\lpa >= -7= XI CmlVW ( 21 ) 
* m=l 

where 

Ci = — sin^s; (a) 

c 2 = -sin^; (6) 

c 3 = cos6»b; (c) 

c 4 = cos 0a (d) (22) 

In Eq.(22), 0oa and #ob are replaced by 0a and 9 B respectively If \ip a > in Eq.(19) is the full system state at time 
t=0, then its evolution is given by 

I 4 

\1> a (t) >= -= V c fe exp(-^ fe V jh)\i) k> (23) 
^ 2 fe=i 

Eq.(23) gives the general case of time evolution. We can study it particularly for our proposed system as a special 
case and discuss the pattern of superposition of wave functions. 

B. Special Case 

To analyze the time development, first consider the special case of n A = n B = k, E a = E B = E at0 m, and 
u> A = lu b = ui. Then, e A = e B = e . It also should be noted that then the two photon modes have the same frequency. 
Since we have assumed that the two photon modes are distinct, the two polarization directions of the mode must be 
perpendicular. Also note that we are not necessarily at resonance, i.e., e is not nessarily equal to zero. 

Continuing, for this special case, we also have q A = q B = q with 



and 6 A = B =9 where 9 is given by 



A 2 + (^) (24) 



sin * = tinV^W (25) 



The time evolution of state of the full system is then 
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\ip a {t) >= -^{-sintfcxpH^t/ A) [IV'i > +1^2 >] +cos6»exp(-^3t/ /i) [\ip 3 > + |V> 4 >]} (26) 

To obtain Eq.(26), we have made use of the fact that Ei = E\ and E4 — E 3l for this special case. The values of E\ 
and E 3 are 



El — E a tom 
E 3 E a tom 

We have not yet assumed the resonant behavior, when the detuning parameter is taken to be zero. 



(27) 



C. Resonant Subcase (€= 0) 

The evolution of \tp a (t) > is particularly simple for the subcase in which the detuning parameter e= 0. Then, Eq. 
(25) yields 



cos 9 = —= ; 
V2 



(28) 



Further, from Eq.(24), 



smtf 



I A 

V2W\ 



(29) 



9=|A|; 



(30) 



E\ — E a tom [1 — 



(31) 



Then, Eq.(26) simplifies to the result 



E 3 = E atom [1 + |A|] 



(32) 



\Mt) >=^exp(-iw/ h) |-A e ^ m |Ai^[|^ > + |^ 2 >j + e -iE atom \x\t/ V[ ^ 3 > + |^ (33) 

To interpret the time changing nature of \ipa(t) >, we expand the two square brackets in Equation (30) in the basis 
>, k=l, 2, 3, 4. So, we can then write 



(34) 



k=i 



with 



h = h = - 



1W\ 



/3 - h — 



V2~ 



..(a); 
.(b) 



(35) 
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Also 



|V>3 > +1^4 >= X^ fc l $fe > 



fe=l 



with 



(36) 



Si = 52 



1 

1 A 

33 = 54 = ^]A| ■ 
Under these conditions, Eq. (33) can be simplified to be 



...(a) 
■(b) 



(37) 



(jgl > + |$2 >) 

A |$ 3 > +|$ 4 > r 



e iE atom \\\t/lJ + e -iB otom |A|t/^' 



|A| 2V2 
which can be represented in angular form as 



e ~Ea t om|A|t/# _ e -iE atom \X\t/H 



(38) 



e -iE atom tm f / E 
\Wa{t) >= ~ < COS 



V2 



h (1$! > +|$ 2 » - zsin ( ^ |A|t ) A (|*3 > +l$4 » > 18!)) 



Eq.(39) shows entangled states similar to those in Eq.(7) in PC. However, our result is more general as we have 
included the non-interacting part of the energy in our model. Moreover, our results could further be generalized to 
n-particle states also. However, Eq.(33) can still be further generalized for the case of nonzero detuning. 

D. Non-Resonant Subcase (£j^ 0) 

In general, the energies of the photon mode will not exactly match the energy difference between the atomic ground 
state and the atomic excited state. That is, there is some detuning and 0.. 

Then it is straightforward to extend the results of Part C to allow detuning. The expensions shown in Eqs.(34) and 
(36) remain valid, but the coefficients f& and gk can be easily generalized to the results 



and 



fi 
h 



9i 

93 



h 
h 



— suit/ 
cos#; 



52 = cos 0; 
g 4 = sin 6. 



(a) 
(b) 



(a) 
(b) 



(40) 



(41) 



Finally, Eq.(39) is replaced by the more general result 



>= z~ tEat ° m t/h {F(q, 6, t)\i> a > +G(q, 9, i)|^>} . 
In Eq.(42), we have employed the definitions 



(42) 
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Also note that 



, 3 

E atom = -Eatom(l + ^ G ); ( a ) 

F(q,9,t) = cos(qE atom t/ /i) - i(cos29)sm(qE atom t/ fi); (b) 

G(q,9,t) = -i(sm26)sm(qE atom t//i). (c) (43) 



|Va > = |^a(0) >= ^ {1*1 > +l $ 2 >> S («) 

l# > --^={|<I> 3 >+|$4>}. (6) (44) 

In Section IV that follows, we will discuss the new effects present for entanglement exhibited by Eq.(42), as 
contrasted with the resonant case result of Eq.(39) 



IV. COMPARISON OF PC RESULTS WITH JC MODEL RESULTS 

To compare the result of Eq.(39) with the results of PC, we compare their notation with ours with the following 
correspondance 

10 > |->|->=|G> (a) 

"' = + ( b ) 

(a+|0>)|- > |->=-L(|$ 1>+ |$ 2 >) ( C ) 

|0 > |+ > |- >= |*3 > (d) 

|0 > |- > |+ >= |$ 4 > (e) (45) 

In the study of PC, the authors set 

E a tom 

A — 1 — 7 

A 

Comparison shows agreement of our results with those of PC (for positive A) with their equation (7) except 

(a) The oscillating factor e -^*°™*/^ f Eq.(39) is missing in Eq.(7) of PC. This is because they are assuming that 
the complete Hamiltonian is the interaction term coupling the atom to the photon modes. 

(b) Eq.(7) of PC does not contain the factor -i in the second term of our Eq.(39). However our results include PC's 
results as a special case. 

We can bring agreement between the second term of Eq.(39) and the second term of Eq.(7) of PC by the following: 
Replace the excited state atomic wave functions for both atoms in the work of PC by — (PC). Here, <&j+ (PC) 
denotes the excited atomic wave functions as used by PC. (This replacement is simply a multiplication by a phase 
factor and hence is an equally valid representation of the excited states). It then quickly follows that the second terms 
of our Eq.(39) and Eq.(7) of PC are identical. 



V. DISCUSSIONS 



We have studied entanglement in the rotating wave approximation in an extension of the model used by PC. 

First, consider the resonant case result (e= 0), disployed by Eq.(39). At time t=0, only the first term of Eq.(39) 
is non-zero. At this time, both atoms are in their ground state (see Eq. (18)). However, at the later time 
t=to=^ ( \\ \ e t — )i on ly the second term is present. Then from Eqs. (17) and (33), the following entanglement 
is deduced: If atom A is in ground state, then we know, with certainty, that atom B is excited, and vice versa. In 
other words, it is impossible for both atoms to be in their ground states, at time t=t . 
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There is a periodic increase and decrease of the strength of the entanglement, as expressed by Eq.(39). The period 
T for a full cycle is T=( |A ff om ). 

Now turn to the more general non-resonant case, expressed by Eq.(42), (e^ 0). Then the statements in the 
preceding paragraphs are modified by noting that maximum entanglement now occurs at time t=t = ( 2q E^ t — )■ 

Now, from Eq.(24), |A|<q, for 0. Hence, t < to, when G^O, and maximum entanglement occurs more quickly 
than in the resonant case. The period of oscillation is now T = ( ) ■ 

To summarize, in this paper, we have studied entanglement of two atoms and two photon states and its time 
evolution. However, we have entered into a model that can be extended to a larger collection of atoms in the presence 
of a larger number of photon modes. Moreover, using the modified form of our master equation, we can calculate 
the dissipation of any energy mode from a cavity. The master equation gives the major source for the dissipation of 
photon energy. However, it does not account for the contribution to loss due to the interaction of the atoms with the 
cavity. This dissipative dynamics of cavity can be derived from the leakage of cavity photons due to the imperfect 
reflectivity of the cavity mirrors. It is usually considered in the JC model that the presence of atoms in a cavity may 
not significantly affect [7] the cavity losses. However, due to the possibility of entanglement, it may not be true any 
more. 
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